Light-driven liquid crystalline nonlinear oscillator under optical periodic forcing 
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An all-optically driven strategy to govern a liquid crystalline collective molecular nonlinear os- 
cillator is discussed. It does not require external feedbacks of any kind while the oscillator and a 
time-depending perturbation both are sustained by incident light. Various dynamical regimes such 
as frequency-locked, quasiperiodic, forced and chaotic are observed in agreement with a theoretical 
approach developed in the limit of the plane wave approximation. 
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The response of nonlinear oscillators to time- 
dependent external control parameter is well-known and 
has been investigated in a wide range of systems includ- 
ing physical [T], chemical Q and biological ones [§]. Pe- 
riodic forcing may lead to entrainment or quasiperiod- 
icity according to whether the ratio fp/ In between the 
forcing frequency (/f) and the natural frequency of the 
autonomous system (/at) is rational or irrational. The 
natural limit cycle behavior of the oscillator may also be 
driven into more complex dynamics. Recent applications 
of periodic forcing such as the control of a chaotic chem- 
ical reaction [3] or the control of a microfluidic droplet 
emitter Q have emphasized the use of this technique to 
engineer controllable systems. 

An optical forcing scheme to achieve a non-contact 
control of a nonlinear oscillatory system is limited to a 
few number of systems. A well-known example is the 
light-sensitive Belousov-Zhabotinsky reaction, which is a 
spatially extended system where a homogeneous time- 
periodic optical forcing allows to control pattern for- 
mation 0. A liquid crystal light valve, where a two- 
dimensional external feedback and an additional qua- 
sistatic electric field have been used, is another exam- 
ple In contrast, a local optical periodic forcing on 
uniformly oscillatingcatalytic surface reaction has been 
investigated in Ref. [8], where the heating provided by a 
focused laser beam acts as the local external perturba- 
tion. To our knowledge the all-optical situation, where 
both the nonlinear oscillator and the time-dependent per- 
turbation are supplied by light only, without need of an 
external feedback, has not been reported yet. 

We propose an all-optical periodic forcing strategy to 
control a liquid crystalline molecular oscillator driven 
by light based on optical orientational nonlinearities of 
mesophases Q . There is no need of external feedback ow- 
ing to the strong light-matter coupling that occurs during 
the propagation of light through the optically anisotropic 
liquid crystal (LC) material. Indeed it is known that a 
LC under an intense laser field can be viewed as a collec- 
tive molecular oscillator whose dynamics depends on the 
light-matter interaction geometry. Very rich director dy- 
namics has been observed and confirmed theoretically at 
fixed light intensity. Among most studied geometries one 



can mention a circularly polarized beam at normal inci- 
dence [lol . , an elliptically polarized beam at normal 
incidence jl2l . IT3| . an ordinary linearly polarized beam 
at oblique incidence ^4*1 or a linearly polarized beam 
at normal incidence having an elliptic intensity profile 
p^ . More specifically, laser- induced nonlinear reorien- 
tation dynamics in LCs has retained some attention in 
the context of transition to chaos [l^, [13, [3 ^^^"^ chaotic 
rotations [l^ [23| . In contrast to all previous studies per- 
formed at fixed intensity, the aim of the present work 
is to explore both theoretically and experimentally two 
representative light-LC interaction geometries when the 
intensity is periodically modulated. 

The observations are performed using the experimental 
set-up shown in Fig. [T] where the linear geometry refers to 
a linearly polarized light impinging at oblique incidence 
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FIG. 1: Experimental set-up. Multi-line Ar"*": Argon-ion 
laser; DP: dispersion prism; M: mirror; A/2: half-wave plate; 
P: polarizer; R-A/2: rotating half-wave plate; BS: beam split- 
ter; A/4: quarter-wave plate; L: lens; NLC: film of nematic 
LC; d: diaphragm; IF: interferential filter for 514.5 nm; D^: 
photodiodes. (a,b) linear and circular geometries; ko and E: 
incident wavevector and electric field. 
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(typically few degrees) onto nematic film with the elec- 
tric field perpendicular to the incidence plane [Fig.[T]Ja)]. 
The circular geometry corresponds to a circularly polar- 
ized light with normal incidence [Fig.[TJb)]. In both cases 
the autonomous system has a limit cycle behavior at the 
onset of a secondary Hopf bifurcation, where the reori- 
entation amplitude |nj^| (nj^ = n — no and no is the 
unperturbed director, which is the unit vector that repre- 
sents the averaged local molecular orientation, see Fig. [J) 
acquires an oscillatory behavior with natural frequency 
/jv- The periodic forcing is then achieved by sinusoidal 
perturbation of the intensity, p{t) — po + Sp sin(27r/i?t). 
Here, p is the intensity normalized to the reorienta- 
tion threshold, which corresponds to the so-called optical 
Freedericksz transition, and Sp is the forcing amplitude. 

We used in the experiment commercial nematic 
E7 from Merck which was placed between two glass 
substrates that were chemically treated to ensure 
homeotropic anchoring (molecules are perpendicular to 
the film walls). The incident beam is focused onto a ne- 
matic film of L = 75 fim thickness using a 150 mm focal 
length lens. The aforementioned general form of p{t) was 
obtained by combining two laser lines (Ai = 514.5 nm 
and A2 = 488 nm) selected from a multi-line linearly po- 
larized Argon-ion laser using a dispersion prism placed 
at the laser output. The total intensity of each of the 
two beams is controlled using independent combinations 
of a A/2 plate followed by a polarizer. The linear po- 
larization of the A2-beam is continuously rotated owing 
to a rotating A/2 plate controlled by an electrical motor. 
The polarizer placed after the rotating A/2 plate ensures 
that A2-beam is polarized along the y direction with the 
intensity p2{t) oc sm^{'!Tfpt) whereas the intensity pi of 
the y-polarized Ai beam is kept fixed. Both beams then 
recombine through a beam splitter to generate an excita- 
tion light field whose total intensity p{t) = pi+p2 (t) is of 
the required form. Note that the forcing frequency (fp) 
and amplitude (Sp) can be adjusted independently. Fi- 
nally, a A/4 plate placed before the lens allows to switch 
between two geometries. Its optical axis is set along (at 
45° of) the y-axis in the linear (circular) case. 

At first the total intensity, which is monitored by the 
photodiode Utot (see Fig.[T|), is increased smoothly from 
zero, setting Sp = 0. The homeotropic state remains sta- 
ble below the Freedericksz threshold. Above threshold 
the system settles either to a stationary distorted state 
(fixed point) in the linear case or to a state of uniform 
precession of the director around the z-axis with the fre- 
quency fp (limit cycle) in the circular case as shown in 
Figs.^bjb'). A further increase of the intensity po leads 
to a secondary supercritical Hopf bifurcation in both 
cases. At the onset of this instability, a new frequency 
/at associated with the oscillation of the reorientation 
amplitude appears as sketched in Fig. [2Ic,c'). Forcing 
experiments are performed slightly above the secondary 
threshold (typically a few percent) ensuring that a sys- 
tem never reaches higher instabilities for both linear UM 
and circular geometry [11| |. 



Both the natural limit cycle ai Sp = and the re- 
orientation amplitude dynamics ai Sp ^ are moni- 
tored by the time-dependent total intensity of the cen- 
tral part of the beam that emerges from the sample, 



where the contribution /; 



4'°*(t) = /r^\i) 

depends on pi and |nj_|. Indeed Ic^''' cx pi at fixed |n^|. 
In addition, the larger |nj^| is, the stronger self-focusing 
effects are and Ic is small when pi is fixed. Therefore the 
Fourier spectrum of Ic^^\t) can safely be associated with 
the one of |n_L|(i) since pi is constant. This is achieved 
by using an interferential filter operating at Ai to get rid 
of the A2 contribution where p2 (t) is modulated with fre- 
quency fp. The photodiode Dc (see Fig.[T]) thus collects 

the signal Ic^^^ that will be further denoted as Ic- We no- 
tice that we were able to fully characterize the director 
dynamics (i.e. its polar and azimuthal degrees of free- 
dom) using a polarimetric analysis of the output beam, 
which is not shown in Fig. [T] However, only the polar 
degree of freedom is of interest in the presented results. 

A rigorous description of the dynamical properties of 
nematics is well-established in a full hydrodynamic ap- 
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FIG. 2: Maps of the dynamical regimes near 2:1 entrainment 
region in linear and circular geometry. (a,a') Lp/,: frequency- 
locked (/ = /f/p); F: forced (/ = /f); QP: quasiperiodic; 
PD: region of period doubling [/ = fp/i^n), n — 2,3...]; C: 
chaos. The director dynamics on the unit sphere is sketched 
below (b,b') and above (c,c') the secondary Hopf bifurcation. 
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proach [24]. The basic equations are those of the di- 
rector and the velocity field, which are coupled with 
the Maxwell's equations that governs the propagation of 
light. The director equations are derived using the torque 
balance condition among the elastic, electric and viscous 
torques. The flow caused by the director reorientation 
was shown to lead to quantitative rather than to qualita- 
tive changes in both geometries p^, ^23] within the range 
of intensity that has been explored here. Thus, the veloc- 
ity is neglected in our study. In the calculations we used 
the plane wave approximation and assumed that the di- 
rector depends only on z and t. Taking into account that 
A -C L the Maxwell's equations were solved under the 
geometrical optics approximation. The latter are conse- 
quently reduced to a set of two ordinary differential equa- 
tions for the amplitudes of the ordinary (o) and extraor- 
dinary (e) waves. Despite the simplification described 
the resulting set of equations is rather cumbersome and 
is not presented here explicitly. The complete bifurcation 
scenario for the autonomous system {ip = 0) is available 
for both linear [2l| and circular geometry [ll| under the 
aforementioned approximations. The models developed 
there might straightforwardly be extended to the case of 
time-dependent intensity p{t). As in previous studies, 
we used different representations to describe the director 
depending on the geometry. Such a choice is dictated 
by symmetry considerations. Then an expansion of the 
director components is performed in terms of orthogonal 
functions which satisfy the boundary conditions n^^yiz — 
0, L; t) = 0. Namely we took n = (sin 6, cos 9 sin ip, 
COS0COS0) with = sin(n7rz/L) and 

6 = ^^^.j^ 0„(i) sin(n7rz/L) for linear geometry, 
whereas n — (sinOcos$,sin0sin$,cos0) with 
= En=i sm{mTz/L) and $ = $o(i) + 

J2'^=i ^n{t) sin[(n -|- l)7Tz/L]/ sm{Trz/L) for the circular 
one. The Galerkin procedure was then used to obtain 
a set of nonlinear ordinary differential equations for the 
mode amplitudes ((/>„, 6m) or $m) which was solved 
numerically using standard Runge-Kutta method. In fact 
it is enough to retain only a few number of modes for 
the director expansion to obtain a good accuracy for the 
calculated director components (better than 1%). Note 
that the Maxwell's equations were solved at each step of 
numerical integration for time t. We also introduce the 

total phase delay A{t) = 2-k/L j^[ne{z,t) — no]dz be- 
tween the o- and e-waves across the whole film (ne,o are 
the refractive indices), which is a global measure of the 
amplitude of reorientation The calculated A(t) will 
thus be compared with the measured Ic{t). 

The calculated maps of dynamical regimes are shown 
in Fig. [2] in the {Sp/ po, Jf/ In) plane, where the same 
material parameters as in Ref. 11] have been used. In 
addition, po is taken approximately 1% above the sec- 
ondary Hopf bifurcation threshold. For the sake of illus- 
tration the results are presented for a typical layer thick- 
ness L — 100 fim and a typical incidence angle a = 5° (for 
the linear geometry only, see Fig.[T]). These maps aim to 
show the generic dynamical behavior of the system and 




FIG. 3: Demonstration of QP, L2/1 and L3/2 regimes in linear 
geometry for Sp/po — 0.05 and a = 3.5°. Gray (black) curves: 
experimental (calculated) power spectra of Ic (A), (a) right 
part of the L2/1 tongue at fp/fN = 2.37; (b) inside the L2/1 
tongue at fp/fN = 2; (c) left part of the L2/1 tongue at 
fp/fN = 1.62; (d) inside the L3/2 tongue at fp/fN = 1.5. 
For QP states (a) and (c) the frequency peaks (1-5) refer to 
(/iv, fp - fN, 2fN, fp, fN + fp) respectively. 

Linear geometry 




FIG. 4: Director trajectories in the {cj>i,(j)2,6i) [(9i, 62, $1)] 
phase space for linear [circular] geometry illustrating the tran- 
sition to chaos at Sp/po — 0.3 [Sp/po = 0.27] for the same pa- 
rameters as for the maps shown in Fig. [2] (a,a'): fp/fN ~ 1.5 
[1.767] (Umit cycle in the L2/1 regime); (b,b'): fp/fN = 1.35 
[1.742] (limit cycle after the first period doubling); (c,c'): 
fp/fN ~ 1.27 [1.691] (hmit cycle after the second period dou- 
bhng); (d,d'): /f//]v = 1.18 [1.632] (chaotic attractor). 



further quantitative comparison with experiments is done 
using the actual experimental values. We note however 
that the conclusions are unchanged under (i) cell thick- 
ness changes as long as the geometrical optics approxi- 
mation is satisfied and (ii) incidence angle changes in the 
typical range 2-8° [for linear geometry] since the observed 
dynamics is mainly dictated by existence of the secondary 
Hopf instability [21]. In both cases we focus our study 
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FIG. 5: Demonstration of the sequence L — > PD^C (cir- 
cular geometry). Power spectra of 7c (experiment) and A 
(calculated) are shown at 5p/po — 0.22. (a,a'): locked state 
at fp/fN =i 2.05 and /f/Zat = 2.06; (b,b'): state within PD 
region at fp/fN — 1.6 and fp/fN = 1.65; (c,c'): chaotic state 
at fp/fN ^ 1.45 and fp/fN = 1.58. 



on the region near 2:1 resonance which is known to be 
strong. An advantage of the system in question is that 
the entrainment region (the so-called Arnold's tongue) 
which emanates from the point (2/1,0) is rather wide, 
which facilitates its observation, see Fig. [2l Experimen- 
tal evidence of entrainment is found by fixing the forcing 
amplitude and scanning the forcing frequency. The re- 
sults for the linear geometry are shown in Fig. [3l where 
the power Fourier spectra of Ic are plotted together with 
the spectra of its theoretical analog, the phase delay A. 
Each panel corresponds to a measurement which lasted 
^ 4000 s, with fp decreasing from the panel (a) to (d). 
The spectra of Ic and A for a state within the 2:1 tongue 
[Fig- [3b)] are characterized by a single frequency, /i^/2, 
as expected. Outside the tongue and for sufficiently small 
dp, the system responds quasiperiodically. In that case 
the spectra of Ic and A are described by means of two 
incommensurable frequencies Jn and fp that generate 
peaks at frequencies ti/at ± mfp, n and m integers, as 
shown in Fig.[31[a,c). The 3:2 tongue is also observed and 
shown in Fig.[31[d). In this case the single frequency fp/3 
determines the dynamics, as expected. For higher forc- 
ing amplitude a forced regime (the F region in Fig. ^ was 
found which has fp as a unique characteristic frequency. 
When starting inside the 2:1 tongue and decreasing the 
forcing frequency (keeping the same value of Sp) , a route 
to chaos via a cascade of period-doubling bifurcation was 
found for moderate to large forcing amplitudes (the PD 
and C regions in Fig. [5]). The resulting typical trajecto- 
ries are shown in Fig. |4]for both linear (at 6p/pQ = 0.3) 
and circular (at Sp/po ~ 0.27) geometries. In Fig. [S] the 



power spectra of Ic (experiment) and A (theory) are com- 
pared for the successive states from the sequence L — > PD 
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FIG. 6: Window of regular dynamics in the chaotic region for 
Sp/po = 0.22 in the circular case, (a) fp/fN = 1.525 [chaotic 
state]; (b) fp/fN = 1.5195 [periodic state with / = /f/3]; (c) 
fp/fN ~ 1.5184 [periodic state with / = /f/6]; (d) fp/fN = 
1.5182 [periodic state with / = /f/12]; (e) fp/fN = 1.5165 
[chaotic state]. 



C for the circular geometry. In particular the tran- 
sition to a chaotic regime is accompanied by an abrupt 
frequency widening around f/fp = 1/2. Such a behav- 
ior is predicted by theory as well [see Figs. [Hb',c')]. In 
addition, the characteristic presence of windows of regu- 
larity inside the chaotic region when fp varies at fixed Sp 
has been found. There, the director dynamics is periodic. 
One of these windows is shown in Fig.[S]for 5p/pQ = 0.22, 
where a /f/3 dynamics suddenly appears when the fre- 
quency is decreased. Period doubling then begins again 
with limit cycles having frequencies /f/6, /f/12, etc. 
and then once again break off to chaos. 

In conclusion, we showed that an all-optical control 
scheme based on periodic forcing of the light inten- 
sity can be implemented in liquid crystalline materials. 
Two light-matter interaction geometries have been used 
for demonstration. Various dynamical regimes such as 
frequency-locked, quasiperiodic, forced or chaotic have 
been shown to result from the optical forcing. The theo- 
retical study is carried out within the general framework 
involving Maxwell's equations and the constitutive ma- 
terial equations. Fairly good qualitative agreement be- 
tween theory and experiments is obtained in the entire 
range of the problem parameters. The polarization sen- 
sitivity of mesophases might also offer the possibility to 
extend present results to a periodic modulation of the 
polarization instead of the intensity. 

This paper is dedicated to the memory of Prof. Lorenz 
Kramer who suddenly passed away on 05/04/2005. The 
preliminary idea of that work was discussed with him. 
We thank Prof. Michael Tribelsky for his valuable com- 
ments. EB thanks the Physics Laboratory of Ecole Nor- 
male Superieure de Lyon in France for the use of its ex- 
perimental facilities. 
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